The quantum Hall plateau transition at order 1/N 
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The localization behavior of noninteracting two-dimensional electrons in a random potential and 
strong magnetic field is of fundamental interest for the physics of the quantum Hall effect. In order 
to understand the emergence of power-law derealization near the discrete extended-state energies 
E n = huj c (n + |), we study a generalization of the disorder-averaged Liouvillian framework for 
the lowest Landau level to N flavors of electron densities (N = I for the physical case). We find 
analytically the large- AT limit and 1/N corrections for all disorder strengths: at N = oo this gives 
an estimate of the critical conductivity, and at order 1/N an estimate of the localization exponent 
v. The localization properties of the analytically tractable N 1 theory seem to be continuously 
connected to those of the exact quantum Hall plateau transition at N = 1. 
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The explanation by Laughlin Q of the integer quan- 
tum Hall effect depends upon an understanding of the 
localization of electrons by disorder in a strong mag- 
netic field. Without disorder, the single-electron eigen- 
states fall into Landau levels at isolated energies E n = 
fiLo c {n + i), n = 0, 1, 2, . . ., separated by the cyclotron 
energy hu c . The effect of a weak disorder potential is 
to displace some weight from the (5-function peaks at E n 
into localized states at nearby energies. Extended states 
persist at energies E n , and the localization length £(E) 
of states at energy E diverges as E — > E n according to 
a power law £(E) oc (E — E n )~ u ■ Experimental results 
on disordered samples j| are consistent with the value 
v w 2.35±0.05 obtained from numerical calculations ||]4|] 
on the lowest Landau level (LLL). 

Current belief is that the quantum Hall plateau tran- 
sition lies in an entirely different universality class from 
the zero-field case, characterized by "two-parameter scal- 
ing" [5|| and a topological term in the tr-model descrip- 
tion [3|. In particular there is now an understanding of 
the minimal features required to obtain numerically scal- 
ing behavior near the transition which is in some 
respects similar to classical percolation, but with differ- 
ent universal properties as a result of quantum tunneling 
and interference . There is also some understanding of 
the apparent insensitivity to interactions of some critical 
indices such as v. || However, relatively little progress 
has been made in finding an analytically tractable de- 
scription of the localization properties near the critical 
energy. The subject of this paper is a generalization of 
the problem to multiple flavors of electron densities, al- 
lowing an analytical approach to the transition. Our dis- 
cussion is based on the Liouvillian approach introduced 
by Sinova, Meden, and Girvin reviewed below. 

The generalization of the disorder-averaged action to 
N flavors of electron densities gives a simple mean- 
field- like theory in the large- N limit. Other large- N 



approaches such as fllO| ] typically generalize the non- 
interacting problem before the disorder average and do 
not obtain v. At first order in the small parameter 1/N, 
we recover anomalous scaling of the localization length, 
i.e., a value for the critical exponent v. Thus we find an 
analytically tractable quantum description whose physics 
seems to connect smoothly in the parameter N to the 
plateau transition (N — 1). The control parameter 1/N 
allows a systematic expansion around N = oo, and the 
large- N limit and 1/N corrections can be found analyti- 
cally for all disorder strengths. 

The localization properties of electrons at energy E 
are contained in the correlation function of the LLL- 
restricted density operators p q : 
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Here Nl is the number of states in the LLL (Nl — > 
oo taken below), £ = y/hc/eB, and (()) indicates the 
quenched disorder average. 

The Liouvillian approach uses the integral over E of 
IL(q,t;E): 
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The key to the approach is that H(q,t) still contains in- 
formation about electron localization but is more eas- 
ily calculated than the fixed-energy quantity H(q,t; E). 
The disorder average for the Fourier transform Tl(q, lo) 
was carried out numerically in and shown to obey 
the form wImL[(g,w) = uj 1 ! 2,11 f(q 2 /oS), where / is an un- 
known scaling function. At mean-field level Q , only dif- 
fusion is found: our model gives a systematic expansion 
beyond this mean-field result. A major goal of this paper 
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is to obtain the prefactor lo 1 I 2v by a 1/N expansion of 

The scaling form for Pi(q,u>) comes about because at 
large R, only states with energies satisfying 
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will have £(E) > R and hence contribute to the correla- 
tion function at long enough times. Here ad is a nonuni- 
versal length set by the disorder potential. Thus the 
integral over energy which gives II(q, lo) is only nonzero 
in a window of size proportional to q 1 ^ . For states delo- 
calized on the scale R= 1/q, ll(q, ui; E) follows ordinary 
diffusive scaling (cjlmll is a function of q 2 /lo). A pref- 
actor to 1 / 2 " rather than q 1 /" gives the scaling function a 
simple form: for q 2 <C lo 



uban.(q,u) « D{Lo){q 2 /uo) 
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with the frequency-dependent diffusion constant D(lo) = 
Douj 1 / 2 ". This form also applies in the classical percola- 
tion limit studied by Gurarie and one of us 

The LLL-projected density operators p q are related 
to the operators r q of the magnetic translation group 
through 113] j5q — e _ 3 £ i r q . The noninteracting LLL- 
projected Hamiltonian is H — 52 q u(— q)f>q, with v the 
Fourier-transformed random potential. Then using the 
commutation relation for the operators r q 



[r q , r r ] = 2i sin 



q A r r q+r , 



(5) 



we obtain the evolution equation for the magnetic trans- 
lation operators: 
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q' 



in terms of the Girvian 
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Then n(g, lo) is a one-body correlation function of the 
Girvian: 



U(q, U ) = ± i {{(q\(u>-Q)- 1 \<i l 



(8) 



Here the states |q) and operator Q are defined through 

(q|0|q'> = G m >. 

Now we take the continuum limit Nl — > oo and as- 
sume a white-noise disorder potential ((v(q)v(— q))) = 
^p-d)(q + q'). The physical frequency loq is replaced 
by the dimensionless combination lo = hu>o£/v. The 
physical propagator is Tl(q, ujq) — —U(q,Lo), with the 
dimensionless propagator U(q, lo) — > — in the clean limit 



lo — > oo. All momenta are dimensionless (scaled by mag- 
netic length i) in the following. 

The result of disorder averaging is an interacting the- 
ory which can be expressed through a functional integral 
over both bosonic (j) and Grassmann ip variables Q : 



U(q,u;) = 
F(u) = 



D(f>D<f> I DijjD^(t) q (j) q e- F ^\ 

rfq(0q'/ , q + V'qV'q) + 

/(lj 2, 3, 4) (/) qi< j6 q2 </> q3 q4 + 

1,2,3,4 L 

2^ qi 0q 2 0q3^q4 + V>qi ^q 2 V'qs V'qi 

The effective interaction from disorder averaging is 
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/(1,2,3,4) = I e -5l*-««*l S «(q a +q 2 - q 3 - q 4 ) 

x sin Qqi A q 4 ^ sin Qq 2 A q 3 ^ . (10) 

The effect of the additional Grassmann variables ( "super- 
symmetry") is simply to eliminate processes with more 
than one density line. We develop a natural generaliza- 
tion of this interacting problem to N flavors of densities 
(N = 1 is the original problem). The N — ► oo limit gives 
an approximate propagator which is similar to the self- 
consistent Born approximation of Q . The utility of the 
large- N approach is that it gives a systematic expansion 
in the parameter 1/N around the diffusive N — > oo re- 
sult, while at the point N — 1 the theory describes the 
exact localization properties of the plateau transition. 

The A^-flavor generalization of the Lagrangian density 
is (only the bosonic part is given, for compactness, and 
the flavor indices i, j run from 1 to N) 



C N = -iu><j>\<t>\ + 
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The real coefficients Cj in (|llj) reproduce the correct 
N — > 1 limit provided that c± + c 2 + c 3 = 1, so there 
is a two-parameter family of generalizations. The vertex 
with coefficient c 3 in (|ll|) does not contribute as N — > oo 
and does not seem to affect scaling qualitatively at order 
1/N, so it is dropped for simplicity. [jl3| We specialize to 
ci = c 2 = 1/2 in what follows: the choice of these coeffi- 
cients equal is "natural" in that the classes of diagrams 
selected by the two vertices have equal weight at order 
1/N with N = 1, as they do in the full theory (i.e., all 
orders in 1/N) at N = 1, For generic c; the theory has 
a C/(l|l) x SO(N) symmetry, which at the point c 3 = 1 
(N decoupled systems) becomes U(N\N). 

In the N — > oo limit, the diagrams with k interaction 
lines which contribute to (4>q4>q) are the diagrams where 
no interaction lines cross, which have the maximum k 
free choices of flavor index (i.e., degeneracy N k ). Only 
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the first interaction term of ([LI]) affects this limit. The 
noncrossing propagator sums these diagrams and satisfies 
an integral equation depicted graphically in Fig. pi 



Il B (q,uj) = 
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The notation IT^g, uj) is that of ||. In the limit g, cj — ► 0, 
q 2 -C w, the noncrossing result IIs(g, w) shows diffusive 
behavior (oTmll = D$q 2 /u>) without the prefactor cj 1 / 2 ". 



FIG. 1. Diagrammatic representation of the large- N prop- 
agator equation (13). The double line is the noncrossing prop- 
agator. 

The zeroth-order diffusion constant Do is a factor of 
V2 smaller here than in the self-consistent Born approx- 
imation of U (the calculation there corresponds to the 
choice ci = 1). The resulting estimate of the critical 
conductivity at the transition, obtained from Do arid the 
exact density of states |14 through the Einstein rela- 

2 

tion, is a xx ~ 0.614-%- compared to the numerical result 

a xx = (0.54 ± 0.04)^ g§ and a xx = ^ 
of Q . Although the large- N estimates of a x 



0.433^ 
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, xx from the 

Liouvillian are of the right order, it should be noted that 
(unlike v) a universal a xx has not yet been found from 
n(g, uj) obtained using numerical diagonalization. 
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FIG. 2. Diagrams contributing to 1/N propagator correc- 
tion. The solid interaction line in (b) indicates one or more 
noncrossing interaction lines. Each diagram of order 1/N has 
one fewer free index choice than the non-crossing diagrams of 
the same order. 

The noncrossing propagator has simple behavior at 
large q, where the argument of the integrand in (|l2|) is 
rapidly oscillatory. In this limit 



tt / v i n B (oo,u) 2 
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or IIb(oo, uj) w — i ■ 



for small i 



The first corrections to the noncrossing propagator 
have k — 1 free choices of index for a diagram of k lines. 
The corrections consist of all maximally crossed diagrams 
of 2 or more lines (using the noncrossing propagator) 
(Fig. |2|a), and in addition possible "rainbows" over the 
maximally crossed portion (Fig. ^|b). The sum of all max- 
imally crossed diagrams can be obtained from the ladder 
sum represented in Fig. |^ since maximally crossed dia- 
grams are related to ladder diagrams after cutting the 
center propagator line and pivoting. 

The sum of all ladder diagrams with incoming mo- 
menta qi , q2 and momentum transfer 1 as labeled in 
Fig. H|, denoted by V(qi , q2? 1), satisfies the integral equa- 
tion 

V(qi,q25 1) = Vb(qi, q 2 ; 1) + 

|^[lUqi+l',q 2 -l';l-l')x 

n B {qt + I', uj)U B {q 2 - I', w)V(qi, qa; 1')] (14) 



where Vo is the original interaction 
Vb(qi,q 2 ;l) = 



qi A 1\ . f q 2 A 1 
sin — - — sin 



(15) 



Note that V has an implicit to dependence. The content 
of (ji~4|), discussed below, is that V has a diffusion pole 
when q 2 ~ qi: 



v(q,q';q' - q) <* 
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q, q^i q, + i 

FIG. 3. Schematic representation of sum of ladder dia- 
grams. The solid block is defined as the sum of one or more 
rungs with the specified total momentum transfer. 



The physics of (16) is similar to that of the weak- 
localization (WL) logarithmic singularity in two dimen- 
sions [0, but with two major differences: the disorder- 
generated effective interaction involves 4 bosonic density 
operators and hence 8 rather than 4 fermionic operators, 
and time-reversal symmetry is broken by the magnetic 
field. The singularity results when the integral equation 
( |l4| ) becomes nearly V = Vo + V, i.e., when the integral 
operator on the right-hand-side has an eigenvalue going 
to 1. Consider a rung of a long ladder diagram. The 
intermediate propagator momenta on the nth rung from 
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an end of the ladder have magnitude proportional to y/n 
(momenta walk randomly in the plane) so most momenta 
in a long ladder can be assumed large. Each rung adds 
two propagators Hb{co,lo) ~ — i +u>/2, and the interac- 
tion averages to —1 + g(g + q') 2 . This gives the estimate 
D x = | in ©. 

The above derivation assumes white-noise disorder in 
the Girvian, although physical white noise disorder cor- 
responds to disorder of correlation length ~ t in Q ||. 
Another assumption is that the lack of an upper momen- 
tum cutoff in the Liouvillian approach does not soften the 
singularity (|l^) to a logarithm. The singularity is found 
numerically to persist without these assumptions, sug- 
gesting that the high-momentum degrees of freedom do 
not change the scaling behavior qualitatively, as in [Ell . 
We call the limit with point disorder and finite cutoff the 
WL limit as these assumptions are standard in that con- 
text. Now we calculate Hmc(<1>u) and obtain a singular 
log to contribution to the diffusion constant. 

The contribution of all maximally crossed diagrams to 
the propagator is 



n B (^) 2 



dq'n s (g>)0/-Vb)(q,q';q'-q) 
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Here the subtraction of Vq removes the first ladder dia- 
gram, which has no crossings and is already in TIb- Inte- 
grating q' in ( |i~6|) leads to a logw contribution in Umc- 
The sum of maximally crossed diagrams (0) already 
shows nondiffusivc scaling (an anomalous prefactor in lu), 
which is modified quantitatively but not qualitatively by 
adding rainbows over the MC diagrams. There is a log- 
arithimic divergence in the sum of maximally crossed di- 
agrams, leading to a diffusion constant 

D(u) = D + ^ + O(N->) 
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using the standard device of re-exponentiation in large- TV 
and RG theories to estimate critical exponents. Hence 



1 

2^ 



+ 0(N- 2 ). 
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The coefficient a is found from the numerical solution 



of (14) and (p~7|) . Fig. || shows sample data from this 
calculation. Numerically D = 0.682, a = 0.18 ± 0.02, 
and therefore v = (1.89 ± 0.1)Y. In the WL limit the 
addition of overloops (Fig. ^b) is found analytically to 
reduce a and increase v by a factor of 2. Preliminary 
numerical results [ p"3| are that the correction is in the 
same direction for the full model. Although this paper 
has focused on the LLL plateau transition, the large-Y 
Liouvillian approach may also be useful for other nonin- 
teracting quantum Hall transitions with a discrete spec- 
trum of extended states. 
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bution 
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4. Numerical results for maximally-crossed contri- 
to diffusion constant D(u), on semilog scale. Here 

-i/n. 
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